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Abstract. The use of carbon nanotubes as optical probes for scanning near-field optical mi- 
croscopy requires an understanding of their near-field response. As a first step in this direction, 
we investigated the lateral resolution of a carbon nanotube tip with respect to an ideal electric 
dipole representing an elementary detected object. A Fredholm integral equation of the first 
kind was formulated for the surface electric current density induced on a single-wall carbon 
nanotube (SWNT) by the electromagnetic field due to an arbitrarily oriented electric dipole lo- 
cated outside the SWNT. The response of the SWNT to the near field of a source electric dipole 
can be classified into two types, because surface-wave propagation occurs with (i) low damp- 
ing at frequencies less than ~ 200-250 THz and (ii) high damping at higher frequencies. The 
interaction between the source electric dipole and the SWNT depends critically on their rela- 
tive location and relative orientation, and shows evidence of the geometrical resonances of the 
SWNT in the low-frequency regime. These resonances disappear when the relaxation time of 
the SWNT is sufficiently low. The far-field radiation intensity is much higher when the source 
electric dipole is placed near an edge of SWNT than at the centroid of the SWNT. The use of an 
SWNT tip in scattering-type scanning near-field optical microscopy can deliver a resolution less 
than ~ 20 nm. Moreover, our study shows that the relative orientation and distance between the 
SWNT and the nanoscale dipole source can be detected. 

Keywords: carbon nanotube, current density, electric dipole, integral equation, near field, scat- 
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1 INTRODUCTION 

Carbon nanotubes, hollow cylindrical rolls of graphene layers, possess remarkable electronic 
properties [1] as well as tremendous mechanical stability and strength [2] that makes them 
attractive for various applications in opto- and nano-electronics [3]. One of the focuses of cur- 
rent research is on the responses of carbon nanotubes to an externally applied electromagnetic 
field [4-6], because carbon nanotubes have prominent absorption resonances in the infrared 
and visible regimes, the resonance frequencies reflecting the carbon nanotube lattice symme- 
try [7]. Moreover, carbon nanotubes are strongly nonlinear optical structures [4,8,9] with an 
ultrafast optical response [10] that could be exploited for ultrafast optical signal-processing 
devices. Furthermore, as research advances, carbon nanotubes are no longer considered as in- 
finitely long. The influence of the finite length through edge effects on the electromagnetic 
responses of carbon nanotubes is being investigated, particularly for carbon nanotubes oper- 
ating as receiving nanoantennas. Reports have been published on isolated single-wall carbon 



nanotubes (SWNTs) [11-13], isolated multiwall carbon nanotubes [14], bundles of carbon nan- 
otubes [15], and carbon nanotube arrays [16, 17]. Actual radio receivers utilizing an electrome- 
chanical modulation of the field emission from a carbon nanotube antenna have been success- 
fully fabricated [18]. 

Most research in this context is devoted to the response of a carbon nanotube to a uniform 
plane wave. But a carbon nanotube's response to the spatially nonuniform near field due to 
a closely located source ought to be substantially different from the planewave response of 
the same carbon nanotube, depending on the relating location and orientation of the carbon 
nanotube and the near-field's source, as is known to be true for other scatterers [19, 20]. This 
can also be inferred from a theoretical analysis of the scattering of a plane wave jointly by a 
metallic nanosphere and an SWNT [21]. More recently, the intensity spectrum of the thermal 
radiation from an SWNT in the near-field zone was found to have additional resonant lines than 
its analog when the SWNT is in the far-field zone [22]. 

An understanding of the responses of carbon nanotubes in the near-field zone is crucial for 
scattering-type scanning near-field optical microscopy (sSNOM) and biomarking. A sharp tip 
placed in the proximity of the sample is used in sSNOM to scatter near fields induced on the 
sample surface by an external focused light beam. The scattered light involves near-field zone 
interaction between the tip and the sample and maps the surface characteristics of the sample in 
terms of a local refractive index and a local absorption coefficient [6]. Thus, a resolution on the 
order of several nm could be obtained by using an SWNT as a tip. Indeed, an resolution of 30 
nm at 633-nm wavelength has already been reported [23]. In the rapidly developing field of the 
terahertz aperture-less near-field spectroscopy and microscopy [24], a resolution of 100 nm has 
been obtained with a tungsten tip [25,26] and 30 nm with a platinum tip [27], and a proposal 
for SWNTs as tips has already been made [13]. 

Another possibility for using carbon nanotubes for SNOM would allow complementary 
spectroscopic characterization of the photoluminescence (PL) properties of nanoscale objects. 
A higher near-field intensity due to the tip changes and amplifies the signal from a source that 
is either a molecule or a man-made nano-object, with detection taking pave in the far-field 
zone [28]. Scanning along the surface produces the PL map with potentially subwavelength 
resolution being controlled by the apex radius of the probe. 

Nanostructures that combine such multiple functionalities as biocompatibility, fluorescent 
signalling, and drug storage and delivery, will advance cancer diagnostics and therapeutics. 
Integration of carbon nanotubes with nanoscale luminescent materials such as quantum dots 
(QDs) appears promising [29]. This integration may be achieved by functionalizing carbon 
nanotubes with either DNA molecules [30] or carboxyl groups [29,31] to which QDs could 
be attached. Integration of carbon nanotubes and QDs substantially affects the luminescence 
properties of QDs, corroborating the experimentally determined energy transfer between carbon 
nanotubes and QDs [31]. 

The foregoing examples amply demonstrate that the near-field responses of carbon nan- 
otubes require a comprehensive investigation. As a first step in that direction, we examine the 
scattering of the near field of an oscillating, point electric dipole in the proximity of an SWNT of 
finite length — to address the effect of PL amplification and to understand the role of the SWNT 
antenna/near-field probe in the formation of the radiation pattern of a point PL source. The out- 
line of the paper is as follows: The boundary-value problem is formulated in Sec. 2, and several 
typical profiles of the current induced on an SWNT are presented and discussed in Sec. 3. Our 
theory is intended to cover the coupling of a broad variety of luminescent nano-objects to an 
SWNT near-field probe. As the energy transfer between a nano-object and an SWNT can be 
either resonant or non-resonant, depending on the excitation frequency, illustrative examples 
are presented in Sec. 3. Section 3.3 is devoted to the special case of an endohedral molecule 
inside the SWNT, the molecule being modeled as a coparallel source electric dipole located on 
the SWNT axis. The resonant coupling of an SWNT placed close to a source electric dipole 



is examined in Sec. 4. Sec. 5 provides a detailed examination of the scattered electric field. 
Examining the contour plots of the far-zone field due to the dipole-SWNT system in Sec. 5.2, 
we determine the spatial resolution that an SWNT tip can deliver in sSNOM. Conclusions are 
provided in Sec. 6. Gaussian units are used, and a time dependence of exp(— iut) is implicit 
with t as time, ui as angular frequency, and i = \f^l. Vectors are denoted in boldface; unit 
vectors are denoted as e x , etc.; and dyadics [32] are double-underlined. 

2 BOUNDARY- VALUE PROBLEM 

Suppose an SWNT of length L is exposed to the electromagnetic field radiated by a current 
density J so (r), where r denotes the position vector. The electric field everywhere must satisfy 
the nonhomogeneous vector Helmholtz equation 

[(V X I) • (V x J) - k% ■ E(r) = ^J s °(r) , (1) 

where k = lj/c is the free-space wavenumber, c is the speed of light in free space, and / denotes 
the identity dyadic. Likewise, the magnetic field everywhere must be a solution of the related 
equation 

[(V x I) ■ (V x £) - k 2 ^ ■ H(r) = —V x J so (r) . (2) 

If the cross-sectional radius of the SWNT is denoted by R cn , the SWNT axis is aligned 
parallel to the z axis of a Cartesian coordinate system (x, y, z), and the centroid of the SWNT 
is designated as the origin of the coordinate system, any point on the surface of the SWNT can 
be specified by 

r cn = R cn (cos 4>e x + sin 4>e y ) + ze z , g [0, 2ir) , z g [— 0.5L, Q.5L] , (3) 

where e x ,y,z are the Cartesian unit vectors. Using the equivalent cylindrical coordinate system 
(p, 4>, z), we have to enforce the satisfaction of the following boundary conditions [33, 34] 

lim {e,, x [H(R cn + 5, <f), z) - U(R cn ~ 5, (f>, z))} 
S— >o 

= —3 eq (z) 1 z e [-0.5L, 0.5L ], 0e[O,27r), (4) 
c 

lim {e p x [H{R cn + S, </>, z) - H(R cn - S, <f>, z)}} 

= 0, z<£[-0.5L,0.5L], 0e[O,27r), (5) 
lim {e,, x [E(i? c „ + 5, (j), z) - E(i? cn - 5, (f>, z)]} 

= 0, z€ (-00,00), 0e[O,27r). (6) 

Here, J eq (z), the surface current density induced on the SWNT's surface S, is a measure of the 
jump in the tangential magnetic field across that surface. 

The surface electric current density is assumed to be independent of 4> and purely axial: 
J eg (z) = J eq (z)e z . For this statement to be valid, the electric field radiated by the current 
density J so (r) needs to be homogeneous along the SWNT circumference. Thus the following 
restriction is imposed on the system under the consideration: kR cn <C 2tt. Based on the spatial 
variations of the electromagnetic field emitted by J so when the scattering SWNT is absent, 
additional restrictions may emerge too. Finally, it must satisfy the edge conditions 

J eq (±0.bL) = 0, (7) 



which express the absence of concentrated charges on the two edges of the SWNT. 



Except in the source region, the electric field can be represented as [32, 35] 



E(r) = E mc (r) +E sca (r) , (8) 

where the incident electric field 



E mc (r) = ^/ d 3 r'G(r,r') • J so (r') (9) 

C l Jyso — 



is due to the source current density that is confined to the region V so lying outside the SWNT, 
and the scattered electric field everywhere is given by 

E sca (r) = ^/ d 2 r'G(r,r') • 3 eq (z'). (10) 
c 2 Js 

In these equations, 

is the dyadic free-space Green function. Equations can be similarly written for the incident and 
the scattered magnetic fields, but are not needed. 

An axial surface conductivity a zz can be obtained to relate the jump in the tangential mag- 
netic field across the surface p — R cn of the SWNT to the axial electric field [33]; thus, 

E z (r cn )=e z -E(r cn ) = ^l, z € (-0.5L, 0.5L) . (12) 

This equation is quite general, except that it does not take into account the effects of spatial 
dispersion as well as excitonic effects; furthemore, we have neglected the contribution of the 
chiral conductivity a Z( f> to the axial current because that contribution is quite small in SWNTs 
[36]. 

Taking the dot product of both sides of Eq. (10) with e z , and making use of Eqs. (8) and 
(12) therein, we obtain 

^1-El nc (r cn )= (JL + k ^Il(z), *e(-0.5L,0.5L), (13) 

where the scalar Hertz potential 

ikJ{z - z 1 ) 2 + 4i?2„ sin 2 (<///2) 



0.5L /.2tt ex P 



U (z) = / dz'ri(z') / 

w J-Q.5L JO 



°- 5L J° yj \z - z') 2 + ARl n sin 2 (</>'/2) 



(14) 



According to our assumptions, El nc (r cn ) changes so slowly along the SWNT circumference 
that it can be assumed to depend on the axial coordinate only; i.e., E l z nc (r cn ) ~ i?™ c (0, 0, z cn ). 
Thereby, Eq. (13) becomes an integrodifferential equation for J eq (z) with the formal solution 



II(z) = C ie ~' lkz + C 2 e lkz + -V 

2ik 



r 0.5L 


f £ ik\z-z'\ 


/ dz' 


/-0.5L 





jeq( 7 ') 



, (15) 



where the constants Ci 2 have to be determined eventually using the edge conditions (7). 
Equating the right sides of Eqs. (14) and (15), we obtain 

0.5L 0.5L 

C 1 e-' lkz +C 2 e* kz + J dz'K,(z-z')J eq (z')^^- J dz' e lk \ z ~ z '\ E l z nc (0, 0, z') , (16) 

-0.5L -0.5i 



with the kernel 



JC(z) 



e ik\z\ 

2ika r .. 



R 



c n 



exp 



ik 



4i?2„ sin 2 (0/2) 



4i? 2 „ sin 2 (</>/2) 



(17) 



Equation (16) is a Fredholm integral equation of the first kind with J eq (z) as the unknown 
function to be determined [37]. It has to be solved numerically for specified E l z nc (0, 0, z), as 
in a predecessor paper [13]. Once J eq (z) has been determined for all z 6 (— 0.5L, 0.5L), the 
scattered electric field can be calculated at any location using Eqs. (7) and (10). 

We chose the source of the near field to be a point electric dipole of moment po located at 
r s outside the SWNT, as shown in Fig. 1; accordingly, 



1-3 



(18) 



(19) 



r- 



J so (r) = -io;po<5(r-r s ). 

Substituting Eq. (18) on the right side of Eq. (9), we obtained 

E mc (r) = fc 2 £(r,r s )-p . 

The electromagnetic field due to an electric dipole contains terms that vary as |r — r, 
r s | -2 , and |r — r 8 | _1 , with the first term dominant in the near-field zone and the third term 
dominant in the far-field zone [32]. There are certain restrictions on po and r s to ensure that 
the z-directed component of the incident electric field varies very little along the circumference 
of the SWNT for any z G (-L/2, L/2). 

In the foregoing equations, we have implicitly assumed that the source incorporates inter- 
action with the SWNT, and is therefore in a steady state in the presence of the SWNT. As 
such, it already includes possible renormalization of the dipole frequency (and lifetime). Such 
quantities can be computed knowing the ones in the absence of the SWNT, but that substantial 
exercise will be taken up in the future. 




Fig. 1 . Schematic for the scattering by an SWNT of the electromagnetic field due to an electric 
dipole of moment po. 



3 TYPICAL PROFILES OF INDUCED SURFACE CURRENT 

Calculations were performed for two types of metallic and one type of semiconductor zigzag 
SWNTs — specified by the dual indexes (14, 0), (15, 0) and (18, 0), as is commonplace [4,7] — 
each of length L = 1 /im. The axial surface conductivity er 2Z was computed using the relaxation 
time r = 3 x 10 -12 s (except for those in Sec. 4) and the overlap integral 70 ~ 2.7 eV, as shown 
elsewhere [33]. 

Two orientations of the electric dipole were considered: parallel to the SWNT axis and 
normal to the SWNT axis. The magnitude p of the electric dipole moment was assumed equal 
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(a) po = 1CP 20 e z esu, r s = 10 e y nm 
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(b) po = 10~ 20 e y esu, r B = 10 e y nm 

Fig. 2. Surface current I(z) induced along the axis of a (15, 0) SWNT of length L = 1 /xm in 
response to a source electric dipole po located at r s in the central transverse plane of the SWNT. 
The electric dipole is oriented either (a) parallel or (b) normal to the axis of the SWNT. 

to 0.01 D = 10~ 20 esu, while the location r s was varied. We illustrate the frequency dependence 
of the electromagnetic coupling between the source and the SWNT in this section by computing 
profiles of the induced current at the following four different frequencies: 

(i) lo/2it = 2.6 THz, the first geometrical-resonance frequency for surface-wave propaga- 
tion, defined by the condition (h 2 — k 2 )L 2 = ir 2 [13], where h is the guide wavenum- 
ber[33, Eq. (58)]; 

(ii) lo/2h = 4.0 THz, an off-resonance frequency; 

(iii) oj/2ir = 5.2 THz, the second geometrical -resonance frequency for surface-wave propa- 
gation, defined by the condition (h 2 — k 2 )L 2 = Ait 2 [13, 33]; and 

(iv) ui/2tt = 500 THz, the frequency of interband transitions between Van Hove singularities 
[39,40]. 

In addition, for the (18, 0) SWNT we chose 

(v) ui/2ir = 1310 THz, a plasmon resonance frequency defined by the condition hw = 2jq 
[13,38]. 



The smallest free-space wavelength A = 2n/k among these five cases is 229 nm. As the cross- 
sectional radius of the (15, 0) SWNT is 0.587 nm and that of the (18, 0) SWNT is 0.705 nm , 
both satisfy the condition kR cn <C 1. 

3.1 Coupling to the near field of a source electric dipole 




z, nm 

Fig. 3. Same as Fig. 2(a), except that the source electric dipole is located at r s along the axis of 
the (15, 0) SWNT (p = 10~ 20 e 2 esu and r s = 510 e z nm). 

The calculated profiles of the electric current I(z) = 2irR cn \ J eq (z)\ induced on the SWNT 
surface by the electric field (19) due to a source electric dipole are presented in 

• Fig. 2(a) for r s = 10 e y nm and po = 10 -20 e 2 esu, 

• Fig. 2(b) for r s = 10 e y nm and po = 10~ 20 e y esu, and 

• Fig. 3 for r s = 510 e z nm and po = 10 -20 e 2 esu. 

The remaining case — r s = 510 e z nm and po = 10~ 20 ej, esu — to complete a quartet turned 
out to be trivial for the following reason: When the electric dipole is situated on the SWNT 
axis and is directed normal to that axis, the z-directed component of the incident electric field 
changes sign along any cross-sectional diameter of the SWNT. Thus, that component of the 
incident electric field cannot be assumed to be invariant across the circumference of the SWNT. 
However, as that component in one half of the SWNT bisected by the meridional plane to which 
the electric dipole is tangential is opposite in sign to that component in the other half, and 
consistently with the fully symmetric one-dimensional approximation allowed by the condition 
kR cn <C 1, we can assume that the coupling between the source electric dipole and the SWNT 
is negligibly small, at least to the first approximation. 

In all other cases (Figs. 2 and 3) the source electric dipole is located 10 nm from the closest 
point on the SWNT, which distance is less than A/20 even at the largest frequency considered. 
Therefore, not surprisingly, the closest part of the SWNT is tightly coupled to the electric dipole 
and sustains a very large surface current density. Apart from that feature, which is common to 
all curves in Figs. 2 and 3, two different regimes of the SWNT's response to the near field can 
be identified. 

The first regime is the low-frequency (terahertz) regime for metallic SWNTs — represented 
by cases (i)-(iii) — wherein the SWNTs support surface-wave propagation with very low damp- 
ing. In this regime the distribution of the surface current density is quite uniform along the 
length of the SWNT, subject, of course, to the edge conditions (7). 

The second regime of the SWNT's response to the near field is represented by cases (iv) and 
(v). These response characteristics are inherent 



• to metallic S WNTs in the high-frequency regime, where the contribution to <j zz from the 
interband electronic transitions becomes essential, and 

• to semiconducting SWNTs in the whole frequency range, as we show next. 

Surface-wave propagation on the SWNT still occurs in the second regime, but with high damp- 
ing [case (v)] and very high damping [case (iv)]. Accordingly, the surface-current-density pro- 
file is strongly peaked [case (v)] and very strongly peaked [case (iv)] closest to the electric 
dipole. 

The surface current profile along the SWNT axis does not depend directly on the free-space 
wavelength A, as the near field of dipole is strongly localized in the vicinity of the dipole. 
In order to confirm this issue, the current induced in the (14,0) SWNT by the electric dipole 
oscillating at the frequency 2.6 THz — which is the frequency of the first geometrical resonance 
in the (15,0) SWNT — is presented in Fig. 4. The dipole orientation and location are the same 
as for Fig. 2(a): po = 10 _20 e 2 esu and r s = 10 e y nm. The current distribution in the 
(14,0) SWNT in Fig. 4 strongly resembles the current distribution in the (15,0) SWNT at 
500 THz — see Fig. 2(a) — though the wavelength is 200 times smaller in the second case. The 
only difference is the maximum value of the induced current I max that is much higher for the 
(15, 0) SWNT at 500 THz than for the (14, 0) SWNT at 2.6 THz. 
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Fig. 4. Same as Fig. 2(a), but for a (14, 0) SWNT. 



The orientation of the electric dipole with respect to the SWNT axis has a significant effect 
on the profile of the induced current. Let us first look at Fig. 2, for the source electric dipole 
located in the central transverse plane of the SWNT. At any point on the SWNT axis, the 
incident electric field has only y- and z-directed components for z G [— 0.5L, 0.5L]. In the 
central part of the SWNT, the incident electric field is primarily axial (z-directed) when the 
electric dipole is parallel, but primarily normal (y-directed) when the electric dipole is normal, 
to the SWNT axis. Therefore, following Eq. (12), at the center of the SWNT the induced surface 
current is not null-valued in Fig. 2(a) but is null-valued in Fig. 2(b). Next, when the electric 
dipole is located on the SWNT axis, as for Fig. 3, the difference is even more prominent as a 
dipole oriented normal to the SWNT axis is not coupled to the SWNT at all. 

3.2 Coupling to the far field of a source electric dipole 

Thus far we have considered the SWNT located in the near-field zone of the source electric 
dipole. It is instructive to consider also the case when the SWNT is located in the far-field zone. 



Although this situation is not helpful to our goal of evaluating the possibility of SNOM applica- 
tions, it gives additional insight on the interplay of the geometrical and frequency dependencies 
of the coupling between the source field and the electronic system of the SWNT. 
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Fig. 5. Same as Fig. 2(a), except for r s = 100 e y /jm. 



The profiles of the induced surface current computed for r s = 100 e y fim and po = 
10~ 20 e z esu are presented in Fig. 5. Two response regimes are evident in this figure, just as in 
Fig. 2(a): one regime covers cases (i)-(iii), and the other spans cases (iv) and (v). The magni- 
tude of the induced current is much smaller when r s = 100 e y than when r s = 10 e y nm, 
because the incident electric field is much weaker in magnitude at any point on the SWNT when 
the SWNT is irradiated by the far field. The profile of the induced current tends to be more uni- 
form for r s = 100 e y fim than for r s = 10 e y nm, because every point on the SWNT lies in the 
far- field zone of the source electric dipole in the former scenario. 

The oscillatory profile of the induced surface current in Fig. 5 for case (v) is due to the fact 
that the wavelength of the surface wave is ~ 200 nm. In order to explain this issue, let us recall 
that a surface wave with exp(ihz) dependence on z can travel on the SWNT, per [33, Eq. (57)]. 
The value of h works out equal to h = 3.1 x 10 5 + i6.4 x 10 4 cm -1 , by virtue of [33, Eq. (58)], 
and the surface wave therefore has wavelength 27r/Re(/i) ps 200 nm. 

3.3 Coupling to a source electric dipole inside the SWNT 

The formulation presented in Sec. 3 can be easily modified to to accommodate the location of 
the source electric dipole inside the SWNT, so long as the z-directed component of the incident 
electric field depends on the axial coordinate only. This situation mimics the special case of 
an endohedral molecule inside the SWNT, the molecule being modeled as a coparallel source 
electric dipole located on the SWNT axis. 

Accordingly, after setting po || e z , r s = z s e z , and z s € (— 0.5L, 0.5L), we have to modify 
Eq. (15) to 



n(z) 



Cie 



-ikz 



C 2 e ikz + — 



1 

2ik 



0.5L 



dz , ) e ik\z-z'\ 



0.5L 



J e <V) 



E l z nc (R cn cos </>, R cn sin <f>, z') 



(20) 



with the tacit understanding that E z nc (R cn cos <j>, R cn sin <f), z') is independent of < 
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Fig. 6. Same as Fig. 2(a), except for the source electric dipole placed at the centroid of the 
(15,0) SWNT (r s = nm). 

The profiles of the surface current induced in the SWNT by a source electric dipole placed 
at the centroid of the SWNT (r s = nm) are presented in Fig. 6. These profiles resemble the 
ones in Fig. 2(a) for the electric dipole placed at r s = 10 e y nm. However, we observe much 
more prominent rise in the induced surface current in the parts of the SWNT close to the source 
electric dipole in the former case due to the higher strength of the incident electric field at the 
SWNT surface. 

4 RESONANT COUPLING IN THE NEAR-FIELD ZONE 

The edge conditions affect not only the profile but also the magnitude of the induced surface 
current, the latter being captured by the average induced surface current 



The computed spectra of I avg presented in Fig. 7 contain several resonances at the frequen- 
cies defined by the condition hL as its, where s depends on the location and the orientation of 
the source electric dipole. When the electric dipole is situated near an edge of the SWNT, s is 
an integer; when the electric dipole is located equally distant from both edges of the SWNT, s 
is either even or odd, depending on the orientation of the electric dipole. 

Even though the resonant amplitudes of I avg in Fig. 7 depend on the location and the 
orientation of the source electric dipole, the full width at half maximum (FWHM) of each 
resonance is primarily defined not by the source characteristics but by the attenuation of surface 
waves in the SWNT. This attenuation emerges in our model through the relaxation time r that 
affects the axial surface conductivity a zz of the SWNT. In order to demonstrate the effect of r, 
we present the spectra of I avg for different values of r in Fig. 8. A decrease in r weakens the 
resonances but does not drastically affect the off -resonance values of I avg . The dissipation of 
the surface waves becomes so high for r < 100 fs that the resonances disappear. 

It may seem from Figs. 2 and 3 that the maximum value I max of the induced surface current 
I(z) grows with the frequency and, in general, is much higher at optical frequencies [cases (iv) 
and (v)] than at terahertz frequencies [cases (i)-(iii)]. However, this is not so. The spectrum 
of Imax is presented in Fig. 9(a) for a (15,0) SWNT when a source electric dipole p = 
10~ 20 e z esu is located at r s = 10 nm. The spectrum of the magnitude of the axial surface 
conductivity of the same SWNT is presented in Fig. 9(b), wherein we can identify two different 
spectral regimes: a Drude regime (frequencies < 280 THz) and a regime of interband transitions 
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Fig. 7. Frequency dependence of the average surface current, defined by Eq. (21), induced in a 
(15, 0) SWNT due to a source electric dipole for three different configurations: po = 1CP 20 e z 
esu and r s = 10 e y nm (black line); po = 10 _20 e 2 esu and r s = 510 e z nm (green line); 
Po = 10 -20 e y esu and r s = 10 e y nm (red line). 
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Fig. 8. Frequency dependence of I avg in a (15, 0) SWNT due to a source electric dipole of 
moment po = 10~ 20 e z esu located at r s = 510 e 2 for three different values of the relaxation 
time t. 

(frequency > 280 THz). These two regimes are separated by a strong dip (around 280 THz). 
The spectral characteristics of I max are completely different in these two spectral regimes. In 
the low-frequency regime, I max varies non-monotonically with frequency even though \a zz \ 
does. In particular we observe a number of resonant lines in the spectrum of I max that are the 
geometrical resonances of the surface plasmons — the same as in Fig. 7. In the low-frequency 



regime (< 10 THz) except for exactly at the resonance frequencies, I max is smaller than at the 
optical resonance frequencies. In the high-frequency regime, the frequency dependence of I max 
follows the frequency dependence of \cr zz \; in particular, two poles due to interband transitions 
are clearly seen in both spectra. 
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Fig. 9. (a) Spectrum of the maximum value I max of the surface current I{z) induced in a 
(15, 0) SWNT by an electric dipole of moment po = 1CP 20 e z esu located at r s = 10 e y nm. 
(b) Spectrum of magnitude of the axial surface conductivity \a zz \ of a (15,0) SWNT. The 
resonances due to interband transitions are indicated by arrows. 



To explain the spectral dependence of I max , we should take into account that, according to 
Eq. (13), the induced surface current density in the SWNT arises from the superposition of two 
electric fields; i.e., 



EL 



. d 2 



IL(z) 



z e (-0.5L, 0.5L) 



(22) 



where the second term within the square brackets on the right side is the electric field of a 
surface wave and H(z) is defined by Eq. (14). The nonmonotonic dependence of I max on a zz 
in the low-frequency regime in Fig. 9 is because the electric field of the surface wave depends 
on the axial surface conductivity too. 
Let us recall that 
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Plotted in Fig. 10 are the z-directed components of the incident and the scattered electric 
fields on the surface of a (15,0) SWNT illuminated by an electric dipole of moment po = 
10~ 20 e z esu located at r s = 10 e y nm. At the frequency (= 2.6 THz) of the first geometric 



resonance in the portion of the SWNT close to the source electric dipole, the z-directed com- 
ponents of the incident and the scattered electric fields are of similar magnitude; see Fig. 10(a). 
However, as the real parts of the z-directed components of the two fields differ in sign — see 
Fig. 10(b) — the total electric field on the surface of the SWNT is small, leading to small I max 
despite the high magnitude of the axial surface conductivity. As the frequency rises, the scat- 
tered electric field on the surface of the SWNT rises (especially at the geometrical resonances 
frequencies), as shown in Fig. 10(c), and leads to higher I max despite a lowering of the surface 
axial conductivity. In the high-frequency regime, the magnitude of the axial surface conductiv- 
ity is rather low (except for the interband transitions) and the surface wave is strongly damped, 
the scattered electric field is generally much weaker in strength than the incident electric field, 
as shown in Fig. 10(d); hence, J eq (z) ps <j zz El nc (r cn ) does not depend on the scattered field. 
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Fig. 10. (a) \E' z nc (z)\ and \E s z ca (z)\ at 2.6 THz, (b) Re [Ef c (z)] and Re [E s z c a (z)} at 2.6 THz, 
(c) \E l z nc (z)\ and \E s z ca (z)\ at 60.9 THz, and (d) |£™ c 0)| and \Ef a (z)\ at 400 THz, on the 
surface of a (15, 0) SWNT illuminated by a source electric dipole of moment po = 10~ 20 e z esu 
is located at r s = 10 e y nm. 



5 SCATTERED ELECTRIC FIELD 

5.1 Scattered electric field near the SWNT 

The scattered electric field in the vicinity of the SWNT was calculated using Eq. (10) for all 
five cases, and two locations and two orientations (as appropriate) of the source electric dipole. 
Figures 11 and 12 show |E sca (r)| for r = ze z + 10e y nm, \z\ < 0.6L. 

Substantial enhancement of the maximum value with increasing frequency is not evident in 
Figs. 1 1 and 12. However, those figures do offer evidence of the effect of dipole orientation on 
the scattered electric field. The plots in Figs. 11 and 12 demonstrate that there are two spatial 
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Fig. 11. Magnitude of the scattered electric field |E sca | computed at r = z e z + 10 e y nm, for 
\z\ < 0.6L, when the SWNT and the illumination conditions are the same as for Fig. 2. The 
source electric dipole is oriented either (a) parallel or (b) normal to the axis of the SWNT. 

regions wherein the scattered field is localized: (i) close to the source electric dipole and (ii) 
near the edges of the SWNT. The localization of the field scattered by the SWNT in the region 
close to source electric dipole could be used in SNOM for the excitation of strongly localized 
electric fields caused by the spontaneous decay of an emitter placed in the vicinity of an SWNT. 

5.2 Radiation patterns of the dipole-SWNT system 

Let us now consider the electromagnetic field in the far zone, with direct contribution from the 
source electric dipole and indirect contribution from the scattered field due to the presence of the 
SWNT. For this purpose, we chose a spherical coordinate system (r, 9, <fi) with origin located at 
the centroid of the SWNT. We define the joint radiation pattern 

f(e r )= lim re~ ikr [E SCQ (re r ) + E" lc (re r )l (24) 
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at distances far from the dipole-SWNT system. 
Equation (10) yields 

0.5L 

lim re- lkr E sca (re r ) = , e / 27rJ ^sinfl f e -ikz cos 6 jeq ^ dz ^ (25) 

kr— ¥oo C J 

-0.5L 




-600 



600 



Fig. 12. Same as Fig. 11, except that the irradiation conditions are those for Fig. 3 (po 
10~ 20 e z esu and r s = 510 e z nm). 



which does not depend on (ft, whereas Eq. (19) yields 
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When po || e z , Eq. (26) simplifies to 
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Let us examine the function f (e r ) when the source electric dipole is oriented parallel to the 
SWNT axis and is located on the SWNT axis (i.e., r s — z s e z )\ then, 
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is independent of eft. Displayed in Fig. 13 are polar plots of the normalized joint radiation pattern 
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in any plane to which the z axis is tangential — for z s = and z s = 501 nm, and for the five 
different cases (i)-(v) delineated in Sec. 3. 

In general, the plots of the normalized joint radiation pattern 7(6*) presented in Fig. 13 for 
a (15,0) and (18,0) SWNT do not have the sin 9 form expected of radiation from point electric 
dipoles. The condition for j(9) having the dipolar form — as in Figs. 13(a)-(e) — is kL <C 1, 
which means that the second term on the right side of Eq. (30) can be ignored. But that is not 
enough. Indeed, in case (iii) for the source electric dipole placed near an edge of the SWNT 
edge, Fig. 13(f) does not have the dipolar form even though the condition kL <C 1 obviously 
holds true. Let us study this case in more detail. 

Shown in Fig. 14 the induced surface current density on a (15,0) SWNT at 5.2 THz when 
Po = 10 -20 e z esu and r s = 501 e z nm. Whereas the magnitude of J eq {z) is symmetric with 
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Fig. 13. Normalized joint radiation pattern 7(0) for (left) r s = and (right) r s = 501 e z nm, 
when po = 10 -20 e z esu. From top to bottom: (a, b) case (i); (c, d) case (ii); (e, f) case (iii); 
(g, h) case (iv); and (i, j) case (v). These calculations were made for a (15,0) SWNT (a-h) and 
(18,0) SWNT (i,j). 

respect to the center plane z = 0, the argument is asymmetric. The integral on the right side 
of Eq. (30) is therefore almost null valued for 6 = ir/2 and 7(71-/2) is small. However, for 
8 ir/2 the phase of the integrand changes due to the presence of e - lkzcos6 1 and the joint 
radiation pattern is enhanced. As the right side of Eq. (30) contains terms proportional to cos 9 
and sin 8, radiation lobes centered about | sin 8\ — | cos 8\ = l/v2 appear. 
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Fig. 14. (a) Magnitude and (b) argument of the surface current density induced on a (15,0) 
SWNT by a source electric dipole of moment po = 10 -20 e z esu located at r s = 501 e z nm. 
The frequency is 5.2 THz. 

At sufficiently high frequencies, the condition kL <C 1 does not hold. However, we still 
observe the dipolar radiation pattern in Figs. 13(g,h) for case (iv). When the source electric 
dipole is placed at the centroid of the SWNT, as for Fig. 13(g), the main contribution to the 
normalized joint radiation pattern is from the first term on the right side of Eq. (30). That term 
is manifestly dipolar. Scattering by the SWNT is small in the high-frequency regime, due to 
the strong damping of the surface wave, leading to ohmic absorption rather than reradiation. In 
contrast, when the source electric dipole is placed near one edge of the SWNT, neither of the 
two terms on the right side of Eq. (30) is dominant. But the normalized joint radiation pattern in 
Fig. 13(h) is still dipolar because, as the induced current is confined largely to the edge close to 
the source electric dipole, that edge acts like an additional electric dipole. Accordingly, we con- 
clude that the presence of an SWNT does not lead to PL enhancement at optical frequencies — in 
good agreement with the experimental results where strong luminescence quenching for Gd-Se 
quantum dots and polystyrene spheres in the presence of SWNTs has been observed [41]. 

Figures 13(i,j) indicate that the joint radiation pattern of the dipole-SWNT system is very 
different from that of the dipole alone, in case (v). This happens because the chosen frequency 
is a plasmon resonance frequency. 

The normalized joint radiation patterns presented in Fig. 13 also demonstrate the crucial 
influence of the edges of the SWNT on PL enhancement. At all frequencies the far-field radia- 
tion intensity is much higher when the source electric dipole is placed near an edge of SWNT 
(right column in Fig. 13) than at the centroid of the SWNT (left column in Fig. 13). A strong 
enhancement by a factor ~ 2000 is predicted at the frequency of the first geometrical resonance, 
i.e., case (i). 

5.3 Resolution with an SWNT tip for sSNOM 

As stated in Sec. 1, the formalism presented in this paper can be applied for sSNOM inves- 
tigations of a PL object using an SWNT tip. In order to exemplify that assertion, we fixed a 
(15,0) SWNT of length L = 1 /zm on the z axis with the centroid of the SWNT serving as 
the origin of the coordinate system, varied the location r s = x s e x + y s e y + 505 e, nm of the 
source electric dipole of moment po = 10~ 20 esu radiating at 2.6 THz frequency, and computed 
|f (e x )\/k 2 po as a function of x s and y s for both po || e z and po || e y . The resulting contour 
plots are presented in Fig. 15. 

When the electric dipole is oriented along the SWNT axis, the contour plot of the joint 
radiation pattern in Fig. 15a is a set of concentric circles centered on the SWNT axis. The 
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Fig. 15. Contour plot of \{(e x )\/k 2 po as a function of x s and y s , when a source electric dipole 
radiating at 2.6 THz frequency is located at r s = x s e x + y s e y + 505 e z nm, and a (15,0) 
SWNT of length L = 1 /im is affixed to the z axis with the centroid of the SWNT serving as 
the origin of the coordinate system. The dipole electric moment is (a) po = 10 -20 e z esu; (b) 
Po = 10 -20 e y esu. 

jointly radiated field decreases in strength rapidly: within 20 nm from the SWNT axis, there is 
a drop in magnitude by one order. When the electric dipole is oriented normal the SWNT axis, 
the contour plot in Fig. 15b shows two images arranged symmetrically with respect the plane 
formed by the SWNT axis and the direction normal to both the SWNT axis and the orientation 
of the electric dipole. The two images arise due to the weak coupling of y-oriented dipole with 
the SWNT when the distance between the two is small. 

In order to study the spatial resolution of the SWNT tip as a probe, the calculations for 
Fig. 15 were repeated but with two z-oriented electric dipoles, one placed at r s i = r c + (x a e x + 
y s e y ) /2 and the other at r S 2 = r c — (x s e x + y s e y )/2, where r c = x c e x + y c e y + 505 e 2 nm 
is the radius-vector of the two dipoles system geometrical center, x s and y s are the interdipole 
separations along the x and y-axis respectively. In the contour plots presented in Fig. 16, the 
two electric dipoles cannot be resolved when the inter-dipole separation is 10 nm. However, 
resolution is possible when that separation is 20 nm, leading us to the conclusion that spatial 
resolution of the SWNT-based probe could be between 10 and 20 nm. 

6 CONCLUDING REMARKS 

Following a standard procedure, we formulated a Fredholm integral equation for the surface 
current density induced on a metallic SWNT irradiated by the electromagnetic field of an arbi- 
trary source located outside the SWNT. The integral equation was solved numerically and then 
used to compute the scattered field. Though we chose the source to be a point electric dipole 
for numerical work, the technique can be used for other sources such as electrically small loop 
antennas [19] and even extended sources such as aperture antennas [42]. 

The relative location and orientation of the source electric dipole influence the profile of the 
current induced in the SWNT as well as the scattered electric field in the vicinity of the SWNT. 
This effect can be accentuated by the resonant excitation of a surface wave on the SWNT. We 
proved this by investigating the frequency dependence of the induced current. Strong spatial 
localization of scattered electric field near one or both edges of the SWNT, particularly under 
nonresonant conditions, should promote the adoption of SWNT tips for scanning near-field 
optical microscopy. 

Carbon nanotubes are strongly nonlinear, the nonlinearity becoming substantial when the 
incident power density is on the order of 10 10 W cm -2 or larger [43]. For a gaussian pulse, 
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Fig. 16. Contour plots of |f (e x )\/k 2 po as a function of x c and y c , when two identical electric 
dipoles of moment po = 10~ 20 e z esu radiating at 2.6 THz frequency are located at r s i = 
x c e x + (y c + y s /2)e y + 505e 2 nm and r s2 = x c e x + (y c - y s /2)e y + 505e z nm, and a (15,0) 
SWNT of length L = 1 /im is affixed to the z axis with the centroid of the SWNT serving as the 
origin of the coordinate system. Inter-dipole separation is (a) y s = 10 nm and (b) y s = 20 nm. 



this power density corresponds to a peak electric field on the order of 10 7 V cm" 1 , which is 
three orders of magnitude lower than the atomic field. This value of the electric field strength 
is of the same order of magnitude as the electric field strength at a distance of 1 nm from the 
electric dipole — see Eq. (19) — with dipole moment po = 30 D = 3 x 10~ 17 esu which is the 
typical value of the dipole moment of a quantum dot [44]. Thus, though the power density in 
the near-field region can be quite low, the electric field strength can be high enough such that 
nonlinear effects become significant. Thus, a comprehensive description of the scattering of a 
near field by an SWNT will lead to the development of near-field nonlinear optics of SWNTs. 
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